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In relation to the BSSN formulation of the Einstein equations, we write down the boundary 
conditions that result from the vanishing of the projection of the Einstein tensor normally to a 
timelike hypersurface. Furthermore, by setting up a well-posed system of propagation equations for 
the constraints, we show explicitly that there are three constraints that are incoming at the boundary 
surface and that the boundary equations are linearly related to them. This indicates that such 
boundary conditions play a role in enforcing the propagation of the constraints in the region interior 
to the boundary. Additionally, we examine the related problem for a strongly hyperbolic first-order 
reduction of the BSSN equations and determine the characteristic fields that are prescribed by the 
three boundary conditions, as well as those that are left arbitrary. 
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I. INTRODUCTION 

In a series of papers Q, 0, Q we have carried out a 
study of the relevance of the projection of the Einstein 
tensor G a b perpendicularly to a timclike surface as candi- 
date boundary conditions for the initial-boundary value 
problem of the Einstein equations. Using the Einstein- 
Christoffel formulation Q| we demonstrated that the van- 
ishing of three of the four components of the projection 
G a be b (where e b is the vector normal to the timelike 
boundary) are the three necessary and sufficient bound- 
ary conditions required for the initial value problem, and 
that they are also equivalent to the three necessary and 
sufficient boundary conditions required for the well-posed 
propagation of the vanishing of the constraints. Because 
of the generality of the methods used to demonstrate this 
result, we believe that such will be the case for any for- 
mulation of the 3+1 Einstein equations that is strongly 
hyperbolic. 

The full significance of G a be b = to the more standard 
3+1 formulation referred to as ADM 5] is not as well 
established, essentially due to the lack of hyperbolicity 
of the ADM equations themselves. Yet one can see that 
failure to impose at least one of the four equations leads 
to constraint violations in the solution at later times 

At this time, we are intersted in investigating the rel- 
evance of the Einstein boundary conditions G a be b = 
to the "conformal-traceless" decomposition of the 3+1 
equations due to Baumgarte and Shapiro G], which is 
developed on the basis of earlier work by Shibata and 
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Nakamura Q, and is known as BSSN. We have two rea- 
sons for interest in this formulation. In the first place, 
it is a hybrid between the ADM equations and the more 
recent hyperbolic formulations, because it represents a 
partial reduction to first order of the original ADM equa- 
tions with mixing of the constraints in a manner similar 
to most hyperbolic formulations (see, e.g., ||). As such, 
it represents an opportunity to investigate the transition 
between ADM and its hyperbolic (well-posed) reductions. 

Our other reason for interest stems from the fact that 
numerical simulations performed using codes based on 
the BSSN formulation appear to be able to run for sig- 
nificantly long times, which makes it, currently, a formu- 
lation of widespread choice in numerical relativity (see, 
e-g-, 9]), up to slight modifications (see, e.g., E El). 

In Section[n]we briefly describe the problem of the Ein- 
stein boundary conditions for the ADM equations with 
the purpose of introducing the Einstein boundary condi- 
tions in terms of the fundamental variables of the 3+1 
split. The Einstein boundary conditions in terms of the 
BSSN variables are introduced in Section IIIII whereas 
their relationship to the propagation of the constraints 
is developed in Section llVl In Section we study the 
analogous problem for the case of a well-posed first-order 
reduction of the BSSN equations. Our first-order reduc- 
tion is strongly hyperbolic, having six characteristic fields 
that are incoming at the boundary, and six that are out- 
going at the boundary. Unlike the Einstein-Christoffel 
formulation and other first-order reductions of the BSSN 
equations, it has two subluminal characteristic speeds in 
addition to light speed and rest. We show explicitly how 
the Einstein boundary conditions prescribe only three of 
the incoming characteristic fields, which is entirely con- 
sistent with the results of Section Hvl and coincides with 
our previous results on the Einstein-Christoffel formula- 
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tion. As in our previous work, we restrict to the case of 
vanishing shift vector throughout. We close in Section lVTl 
with some remarks. 

II. THE EINSTEIN BOUNDARY CONDITIONS 

IN TERMS OF THE INTRINSIC AND 
EXTRINSIC CURVATURE FOR THE CASE OF 
VANISHING SHIFT VECTOR 

This section develops a brief review of the boundary 
conditions of the Einstein equations in the ADM formu- 
lation. In our minds, the ADM case acts as the basis for 
most other formulations that are derived from it by linear 
combinations with the constraints and by the addition of 
any number of first-order variables. 

Throughout the article we assume the following form 
for the metric of spacetime in coordinates x a — (x l , i) in 
terms of the three-metric 7^ of the slices at fixed value 
oft: 



surface e b can be used to project the Einstein tensor as 
G a be b in order to obtain the components that have no sec- 
ond derivatives across the boundary. As we have shown 
previously , the vanishing of these can then be imposed 
as conditions on the boundary values for the fundamental 
variables. To fix ideas let's choose a boundary surface at 
a fixed value of x. Thus e b — g bx — (0, j lx ) up to scaling. 
Consequently, G a i,e b = G x , so any linear combination of 
the components of the Einstein tensor with a contravari- 
ant index of value x will be suitable. Explicitly we have: 

G% = -^Y** ((ln 7 ),« -l k %k,l) - KD x a + K x k D k a 

+ a { 1 kl V a K x + 1 lx V lk K k ) (5a) 

G x = K - K% _ 1 + Rij + r2 + K 
a 2 

+R X + - (D j D ja - D x D x a) (5b) 



ds 2 



-a 2 dt 2 + r )ijdx l dx l 



(1) 



where a is the lapse function. The Einstein equations 
G a b = for the four-dimensional metric are equivalently 
expressed in the ADM form pj: 



Jij 



= -2aK, 



(2a) 

a (Rij - 2KuK l s + KKij) - DiDja, (2b) 



with the constraints 



C = - (R - K i} + K 2 ) = 0, 



Ci 



DiKU - DiK = 0. 



(3a) 
(3b) 



Here an overdot denotes a partial derivative with respect 
to the time coordinate (d/dt), indices are raised with the 
inverse metric 7 y , Di is the covariant three-derivative 



consistent with 



7ij 



Ra is the Ricci curvature tensor of 



jij, R its Ricci scalar, Kij is the extrinsic curvature of 
the slice at fixed value of t and K = j^Kij. Expressed 
in terms of the Einstein tensor, the constraints © are 
related to specific components in the coordinates (x l ,t): 



C = a 2 G tt 

d = -OTfijCP*, 



(4a) 
(4b) 



where (|4b|l holds only for vanishing shift vector. The 
constraint character (the absence of second derivatives 
with respect to time) is a consequence of the fact that 
the only components of the Einstein tensor that appear 
in 10} have a contravariant index of value t. In geometric 
terms, @ are linear combinations of G a bn b = where n b 
is the unit normal to the slices of fixed value of t, and is 
therefore given by n a — g ah nb = —ag at = Sf/a. 

Similarly, at any boundary given by a fixed value of 
a spatial coordinate, the normal vector to the boundary 



K 



qx y_ 



KK X + R X - -D x D y a 



K' 
a 



G x = Z -+KK 



R x - —D x D z a 
a 



(5c) 
(5d) 



Here T k lj = (l/2)j kt (j u , 3 + - and the time 

derivative of the components of the extrinsic curvature is 
applied after raising an index, that is: K* = ( , y' lk Kf.j),f 
The reader can verify that R x and R x do not involve sec- 
ond derivatives with respect to x of any of the variables 
and that the combination R x — ^R doesn't either 0. 

Since Eqs. © represent the vanishing of four indepen- 
dent components of the Einstein tensor, they must be ex- 
pressible in terms of linear combinations of the evolution 
equations J2Jl and the constraints 0. For this reason, it 
could mistakenly be thought that they must be satisfied 
identically by the solution of the evolution equations with 
initial data satisfying the constraints. Such would be the 
case only in regions of spacetime where the constraints 
are satisfied, as well as the evolution equations. Since 
the constraints are only imposed on the initial data, it 
is not at all clear that they are satisfied everywhere. A 
brief calculation that involves taking a time derivative of 
the constraints and using the evolution equations yields 



C = ad'C, 
Cj = adiC ^ 



(6a) 
(6b) 



where . . . denote undifferentiated terms. This system 
is strongly hyperbolic [2| because it has real charac- 
teristic speeds and a complete set characteristic fields. 
With respect to the unit vector £ l = 7" /^/"/ xx , normal 
to the boundary, the constraints C v and C z travel with 
zero characteristic speeds, whereas C = C ± C x / Vt 11 
travel with speeds ±a, respectively. Only C y ,C z and 
C + C x I yj-f xx are vanishing at the boundary if they are 
set to zero on the region of the initial slice interior to 



3 



the boundary. Because C — C x / yj-f xx is incoming, it is 
not vanishing at the boundary be virtue of the evolution 
and the initial data alone, therefore it is an example of 
what we refer to, colloquially, as an "un-preserved" con- 
straint, and must be made to vanish at the boundary by 
an appropriate choice of boundary conditions. 

One may not, thus, dismiss the boundary equations ijjjj 
as a matter of principle. On the contrary, it is enlighten- 
ing to see exactly how the boundary equations relate to 
the constraints. 

To start with, by inspection one can clearly see that G x 
and G x are identical to the evolution equations for the 
mixed components Ky and K x of the extrinsic curvature 
(namely: the components j — y and j = z of Eq. (|2b|) 
multiplied by 7 11 ). Thus two of the four boundary equa- 
tions are identically satisfied by the solution of the evo- 
lution equations (irrespective of the initial data) . That is 
not the case with G%. If one uses the evolution equations 
for and for the trace K that follow from (|2b|) to elim- 
inate the time derivatives from G x , one is left with the 
Hamiltonian constraint C. Finally, using the definition of 
the extrinsic curvature 12 a I) in the expression for G x it is 
straightforward to see that G x = aj xl Ci. Summarizing, 
if we represent the evolution equations (|2a[l and (|2b|) in 
the form £?j = and £^ = respectively by transferring 
all the terms from the right into the left-hand side, we 
have 



G? 
G x z 

GX 



aC * + l^^ {dl£ J k _ ^ k d t £] k ) (7a) 



i XJ £f z 



C + -(7 W 5h 
a 



7' 



(7b) 
(7c) 
(7d) 



This means that G x and G x are equivalent to C x and C 
respectively (modulo the evolution), whereas G x and G x 
are equivalent to zero (modulo the evolution). One can 
see that imposing either G x = or G x = is equivalent 
to a boundary condition of the form C~ = BC + with a 
constant B on the system of evolution of the constraints. 
This is a valid boundary condition in the sense that it 
preserves the well-posedness of the system of evolution 
of the constraints 12]. Thus, of the two components 
G x — or G% = 0, one represents a necessary bound- 
ary condition. The other one (or a linear combination 
of G x = and G x = not equivalent to the first one) 
becomes redundant because it represents a condition on 
the outgoing characteristic field C + , which is already de- 
termined at the boundary by the initial values. 

Because the ADM evolution equations are not strongly 
hyperbolic, perhaps nothing else can be said about the 
role of the equations G a bG b — as boundary conditions 
for the ADM equations, except that failure to impose the 
nontrivial one leads to constraint violations with guaran- 
teed certainty. 

In the next Section, we examine the Einstein boundary 



conditions as applied in the formulation of the Einstein 
equations that appears in [||. 



III. BOUNDARY CONDITIONS FOR THE BSSN 
FORMULATION WITH VANISHING SHIFT 

In the case of vanishing shift vector, the formulation 
of the Einstein equations referred to as BSSN consists 
of the following evolution equations 



la 



= -2aA t 

= -SL K 
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A i:j = ae" 4 * 



(8a) 
(8b) 



+2r (iTi} k i + r iTkij — 2DiDj. 
-\l l3 (f k , k +T ku {2T tkl +t m )) 



3 a 



TF 



+KAij - 2A il A l j (8c) 
K = ^DiDja + aiAijAV + \k 2 ) (8d) 

f* = 2a(fVi fcj -^7 ij Xi+6i^,i) 

-2i' y a,j (8e) 
for the 15 variables 



Hi 



7^1n(det7„-) 



K = f j K i: 
4 = e-^iK^-^K 



(9a) 

(9b) 
(9c) 

(9d) 

(9e) 



which mainly separate out the determinant of the three- 
metric and the trace of the extrinsic curvature in order 
to evolve them as fundamental variables in their own 
right. But additionally, three first order variables are 
introduced, T l , corresponding to certain combinations of 
first derivatives of the metric. This has the effect of mod- 
ifying the principal symbol of the system of PDE's. 

For this system, we will express the boundary equa- 
tions iJSJ directly in terms of the fundamental variables, 
in such a way that no second x— derivatives of 7^ or (f> 
occur, nor first x— derivatives of T l , Aij nor K occur. 

We start with the simplest one, that is G x , which con- 
tains second derivatives in the term R x . 

To start with, we have Rij — Rij + Rf^ where Rij is 
the Ricci tensor of 7^ and 

Rt = 



-2D i D j cj)-2Zf ij D l bi<l> 



(10) 
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So R* = i^R y] = e-^r j (Ryi + <) = e-^(R- + 
^ XJ Ryj)- As far as the principal terms are concerned, we 
have 



r j R 



-2D X D„ 



/; _^-r... (11) 



Also as far as principal terms are concerned, the Ricci 
tensor of a metric of unit determinant is 



Ri 



1 



(12) 



We can treat G x with a similar procedure. By analogy 
with the reader should have no trouble to see that 

K x = * ~ 7 7 )lxm,ly 

1 

'^xl-zm z,xm~,zl\~, 



But since 



-i"., • ... 



(20) 



Thus 



: -W 3f r kl (%3M-%l,jk-7jk,yl) + ••• (13) 



2 
1 

2 
1 



o(7 7 -7 7 )7ym,fe; 



~,xm~kl ~ 
+ 2^ 7 7fem,J/i + 



(14) 



The last term can be substituted in terms of the first- 
order variable T l : 



1 ~,xm~kl ~ -pa: , 

7J7 7 7/cm,yi = 2 '^ + ' 



(15) 



Additionally, by expanding the three values of the in- 
dex k in the first term one can see that k — x van- 
ishes identically, so that the first term contains no second 
a;— derivatives of jij as it stands, and we have explicitly 

r>x ^/~xl~yrn ~xm~yl\~ 

R y = 2^7 T "7 T )lym,ly 

+\(T l T m -r m r l )iy m M 

+ir% + ... (16) 



On the other hand, K* = j xm K my = j xm (A my 
\% m K) = A x , so that 



K x — A x 



(17) 



Finally, thus, the principal terms of G x are written out 
in the intended form: 



G'„ 



A* 



/I 

y + e -i<t>l_^xl~ym _ ~ X m~yl^~ 



\2 



ym,ly 



+ ^(7 7 -7 7 hym,lz + 2 1 'V 



-2D X D 



y 



(18) 



By a completely analogous calculation we also have 
A x / 1 

z = ^ e (o(7 7 -7 7 nzm,ly 

+ 2"(7 7 -7 7 )7»,k + 2 r ' z 



(19) 



= f *,* + ... 



(21) 



then the term with an x— derivative of f x in R x cancels 
out with a term in — ^R, so that: 

k-\r = l(r l i ym -r m i yl )%m,iy 



1 ^ (z.xl^^m z,xm~zl\~ 
+ 77(7 7 -7 7 hxmjz 



1 

2 ( 



Additionally, we have K x — A x + ^K, so that 



AI =A% + -K- 



(22) 



(23) 



Collecting the relevant results, the principal terms of G\ 
are finaly written out in the intended form: 



Gx 



A x 2 • /I 

+ K + e 7 ™7 7 )7 



a 3a 
1 



V2 



1 



1 f^,xl~zm ^,xra~zi\~ l T*y 

+ 2(7 7 -7 7 )7a;m : ; z - 2 1 
+2D y D y (j) + 2D Z D Z 



xm,ly 



(24) 



For the remaining boundary equation, G x , we observe 
that 

Gt = -i7 xi ((m(det 7fc 0), it - 7 W 7iw) + ■ • • (25) 

which, with the simple substitution of 7^ = e 4 ^7ij, be- 
comes 

G% ^^e-^(8T l k 1 +r l ,i) + --- (26) 

which we take as the intended form without further sub- 
stitutions. 

In principle, thus, we have four boundary equations 



G x = 0, 

G x — 0, 

Gy = 0, 

G x = 0, 



(27a) 
(27b) 
(27c) 
(27d) 
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with G X ,G X Z ,G X ,G X given by JHH, and 

as candidates for nontrivial boundary conditions for the 
evolution equations JHJ. Notice that, unlike the ADM 
case, none of the four equations is manifestly identical 
to any of the evolution equations (JSJ, the reason being 
that all the evolution equations for the traceless extrinsic 
curvature density Ay contain second derivatives of 7ij 
as part of the Laplacian operator j kl d k di. This makes 
it a nontrivial task to figure out whether the boundary 
equations are identically satisfied by the solution of the 
evolution equations with constrained initial data. In the 
following Section, we argue that at least three of them 
are necessary, on the basis of their role in ensuring the 
propagation of the vanishing values of the constraints. 

IV. RELATION TO CONSTRAINT 
PROPAGATION 

The scalar and vector constraints C and C L of can 
be expressed as follows in terms of the BSSN variables: 

C = e -4 * (f k , k -8D l D l (t)-8D l (j)D l (t) 

+ \l lm 7 jk ,l (7ifc.m-27im.fc)) 

-AijA* - ^K 2 , (28a) 

_ - i _ 2 
C t = y k A k ^-T k A lk --^ jKl A 3k --K,i 

+6i^„ (28b) 

In addition to the scalar and vector constraints C and C,-, 
an additional constraint, 

gi = r+^,j (28c) 

must be imposed on the initial values, i.e., 

Q l = (29) 

in order for the BSSN equations JHJl to yield a solution 
of the ADM equations J2J- Whether the constraints C,Cj 
and Q l are satisfied by the solution of the evolution equa- 
tions © is determined by the equations for C,Ci and Q l 
that are implied by © . The easiest propagation equation 
to obtain is that for Q % , because Q % is a trivial constraint, 
whose propagation should be exactly zero by construc- 
tion unless some constraint had been used in establish- 
ing the evolution equations for F* or for 7^. Since the 
vector constraint is necessary to write Eq. l|5c|) in such a 
form, the evolution of Q l naturally depends on the vector 
constraint, as follows: 

Q l = 2aj ij C j (30) 

This equation has no principal terms for the purposes of 
the structure of characteristics, so we can write it as 

Q l = ■ ■ ■ (31) 



where . . . represents undifferentiated terms. Addition- 
ally, since the time-derivatives of the variables that ap- 
pear in the expression (|28a|l for C will not contribute 
any terms that would have derivatives of the constraints 
themselves, we have 

C = ... (32) 

as well. Finally, the only term that contributes deriva- 
tives of constraints to Ci is ^ k A kiy j, yielding 

d = ^Y~ *fikD l D l Q k + ^C, i + ... (33) 

As it stands, the system of evolution equations (f 3 II) . 132|) 
and (|33() is not manifestly well posed, not even mani- 
festly hyperbolic, because of the presence of second space 
derivatives in the right-hand side as a direct consequence 
of the addition of the three first-order variables T\ In 
order to facilitate the study of the properties of propaga- 
tion of the constraints, however, we may define a set of 
six first-order variables for the constraints that allow us 
to reduce the differential order of the right-hand side: 

l ik =luG l , k (34) 

These are not to be thought of as additional constraints 
to be imposed on the initial data of the BSSN evolution 
equations. Rather, they will be automatically satisfied 
by any set of initial data that satisfies Q % = 0. But they 
are useful because, by adding them as new variables, the 
propagation of the constraints takes the form: 



C = ... (35a) 

= -C H +^—j kl l m + ... (35b) 
2 

l ik = 2aCi, k + . . . (35c) 

Q l - ... (35d) 



This first-order system of 13 equations is strongly hy- 
perbolic in the sense that it has a complete set of null 
eigenvectors. This means that for each arbitrary covec- 
tor of unit length CiCi7 y = 1 at any given point, there 
are 13 characteristic fields, which travel with real speeds. 
One can easily see that the characteristic speeds of this 
system are 0, and ±a, with multiplicity 7, 3 and 3 respec- 
tively. Thus seven characteristic constraints are "static" 
in the sense that they propagate along the direction of 
evolution d/dt. Of the six non-static characteristic con- 
straints, three are outgoing at the boundary: 

+ Z, = d + l -i k T lk + (36) 

and three are incoming: 

~Zi = C t - X -i k X lk - ~&C (37) 

This indicates that, even if the 13 constraints are set to 
zero on the initial slice, they would not remain zero in 
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the future of the initial slice unless the incoming con- 
straints are set to zero, either directly or indirectly as 
functions of the outgoing constraints. One way to ac- 
complish this would be by means of boundary conditions 
of the form ~Zi = Lij + Zj with constant coefficients L^ , 
which is a necessary condition in order to preserve the 
well-posedness of the propagation of the constraints ^3 • 

Our next goal is to show that the boundary equations 
G x a = with G x a given by JTHJl , EU , (jSJ and JgjJ are 
equivalent, in a sense, to such boundary conditions. This 
would indicate, in effect, that the boundary equations 
G x = are instrumental in preserving the constraints. 

More precisely, we wish to show that the boundary 
equations G x = are related to the non-static constraints 
through linear combinations with the evolution equations 
in the following manner: 



G x 

GX 

X 

G* 
Gi 



-aC x 



(38a) 



}_prnj _ }_„xjcA 
2 5 a l & vi 



-r j ^ + ^£ K (38b) 

a J 6a 



(38c) 
(38d) 



where £$ , £? k , £^ , £ K represent the evolution equations 
(|8b"l) . (|Sa|l . (Scjl and (|8c[j) where all the terms in the right- 
hand side are moved to the left. We can think of the re- 
lationship between the boundary equations and the con- 
straints through the evolution equations as an equiva- 
lence: 



G x ~ 


^ -aC x 


(39a) 






(39b) 


G x ~ 




(39c) 


G x ~ 


- --£ m 2 

2^ -L-zm 


(39d) 



modulo the evolution. 

In order to prove this, we start with the boundary 
equation G% = 0. We have, explicitly from (|5c|) . 



~,xl cA AX n „- 

7 £ x i = A x -ae 



X 



1 



1. 



-2D X D, 



2 - 
3 



D k D 



and also 



k 

(40) 
(41) 



Therefore, adding the evolution equations (with ap- 
propriate factors) to G x as given by Eq. (|24p. we have 



immediately 
-7 



G x + -fitt - ^£ K 



'XJ 



(42) 



-u 



\1 J l lxj,ki + (7 7 - 7 7 )7^m,; a 
1 8 

+ (7 7 ~7 7 HxmM - g- 1 ,k+^-L> -UkW 



-T 3 jk*r\ 3 ) + 



(43) 



The three indicated terms in second derivatives of 7y 
manifestly combine to yield 

(7 7 - 7 7 )7xm,Zt, 

+ (7 7 -7 7 )7xm,/z 

+7 7 7aj,fc2 =7 7 lxi,km- (44) 

Thus (|42|l reads, equivalently, 

1 - 2 - / 
a + -7 - — I = (j (7 7 x ;,fe m - 7 fc:E l , m J 

-i(r fc , fc -8^^))^- 

+ ... (45) 

Both indicated terms are directly expressible in terms of 
constraints now. The second term in parenthesis (with 
the exponential factor included) is manifestly equal to 
the Hamiltonian constraint C, whereas the first paren- 
thesis immediately leads to one of the new first-order 
constraints: 

7 TxZ.fcm — IkxL ,rn = ~ Ikx (f fc + J J ' ,j ) , m + . . . 

->k , 



^IkxQ jrr> 



(46) 



Thus Eq. Ij38b|l is verified. 

We now show how to prove l|38c|l . From l|Sc|) we have 



a 1 

+4D x D y ^j +A X + ... (47) 
Adding this to G x as given by (|18|l we manifestly obtain 

Gx i ^~ xj cA ( /~xl~ym ~xm~yl\~ 

y + -7 J £yj = [h 7-7 7 n y m,i y 

+ (7 7 -7 7 nym,iz 

+T 3 l lm ly 3 ,l m -l X3 lykT\ 3 ) — 

+ ... (48) 

By expanding the index I = x,y,z in the third line of 
the preceeding equation some cancellations take place, 
leading directly to 



1 ~ , 

Gx i „,xj cA {z,xlz-,kmr, x.xj x, "P^ \ i 

„ + -7 J t yj = — — (7 7 7 ym ,fcz - 7 J 7 yfc l 



(49) 
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The terms in parenthesis in the right-hand side are now 
expressible in terms of the new first-order constraint: 

ffw-f'^ri = -r j %k(7 kl ,i+r k ) :j + ... 

= + (50) 

Thus Eq. (|38c|l is verified. A completely analogous pro- 
cedure (with z in the place of the subscript y) leads to 
the verification of Eq. (|38d|l . 

It only remains to verify Eq. (|38a[l . From l|8bjl and H8a(l 
we have trivially 

e+ = 4> + ... (51) 

and 

£? k = Ilk + ■ ■ ■ (52) 

Adding these (in the manner specified by ||3Sa}) to Gf 
as given by (|26|l we have immediately 

G x t +\i xl (&* ,i 

= e- 4 *a (^r l K,i -l xl l kj A kl ,^j +■■■ (53) 

The right-hand side of the preceeding equation is man- 
ifestly expressible in terms of the vector constraint C;, 
with the consequence that Eq. (|38a|l is verified. By 
the set of equations i|38|) . thus, imposing the equations 
G x = on the boundary is equivalent, in a sense, to 
imposing the indicated combinations of constraints on 
the boundary. We now show that such combinations of 
constraints are equivalent to the incoming characteristic 
constraints. 

By Eqs. ® and ^3, using & = (1,0,0)/^/^, 

we have 



G x t - 


- -f r f 


'~3k + + Z t ) 


(54a) 


L ' x r 




- 


(54b) 


SIX 




- + Zy) 


(54c) 


G x , 




- + z z ) 


(54d) 



The four components G x a are equivalent (modulo the evo- 
lution) to linear combinations of the incoming and outgo- 
ing constraints with constant coefficients. Thus imposing 
G x = is equivalent (modulo evolution) to setting well- 
posed boundary conditions on the system of propagation 
of the constraints (~Zi = L\ + Zj). The propagation of 
the constraints requires only three boundary conditions. 
Accordingly, three linearly independent combinations of 
the four equations G x a = are necessary in order to en- 
sure that the constraints will remain satisfied by the so- 
lution of the evolution equations JSJ with constrained 
initial data. A fourth lindearly independent one should 
be satisfied identically. 



This means that by imposing three (linearly indepen- 
dent combinations) of the four equations G x a — on the 
boundary values of the evolution equations (jHJ) the prop- 
agation of the constraints is enforced. 

At this level there is no preferred set of three linear 
combinations out of the four equations G x = 0. This 
freedom is an advantage, because it may easily accomo- 
date requirements that may stem from expectations of 
numerical stability. We do not concern ourselves with 
such considerations in this work. 

The question arises as to whether the fourth equa- 
tion should be used on the boundary values. For the 
sake of argument, we can take the fourth equation to 
be G x + orf m Gf — 0, since this equation is equivalent 
to a linear combination of purely outgoing constraints 
( + 2i^f xl = 0). If one looks at this question from the 
point of view of the constraint-propagation problem 135|l , 
clearly imposing the fourth equation would be redundant, 
given that the outgoing constraints are made to vanish 
at the boundary simply by picking vanishing initial val- 
ues. Yet the fourth equation is not inconsistent with the 
choice of vanishing initial values for the outgoing con- 
straints. 

From the point of view of the evolution equations JSJ, 
however, it is not at all clear that the fourth equation 
should not be imposed. This is because, as it stands, the 
system of evolution equations © is not st rong ly hyper- 
bolic in the standard sense (see, however, [13|])- There- 
fore, there is no clear classification of the 15 fundamen- 
tal variables into static, incoming or outgoing fields. As 
a consequence, to our knowledge, the actual number of 
boundary conditions required by the evolution equations 
(JHJ, being not less than three, remains unknown. 

In the next section we reduce the BSSN equations 
fully to first-order form. The resulting first-order sys- 
tem is strongly hyperbolic in the standard sense and the 
boundary-value-problem can be solved in the standard 
manner. 



V. EINSTEIN BOUNDARY CONDITIONS FOR 
A WELL-POSED FIRST-ORDER REDUCTION 
OF THE BSSN EQUATIONS 

In order for Eqs. iJSJ to reduce to first-order form, the 
spatial derivatives of <p and those of jij not already ac- 
counted for by r l need to be added as new variables (a 
total of 15 new variables). They are 

Qi = (t>,% (55a) 
3 1 

Vijk = Jij,k ~ ^lk(tlj)l,ml lm + ^jlkl,ml lm (55b) 

These new first-order variables have been used before 
in order to find a symmetrizable hyperbolic first-order 
reduction of the BSSN equations by linear combination 
with the constraints ^4|- The difference here is that we 
do not involve combinations with the constraints. By 
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using 7 U to raise the first two indices in i|55b(l we have symbol has thus the form: 



V« k = -f', k +^7 J >\ m -lfijlrn m ~ ki (56) 



One has V^j — automatically, so the set of variables 
V 11 \ does not include the three variables T\ Addition- 
ally, V %3 \ is traceless in the first two indices, also by 
construction (V^klij = 0). Thus V %3 \ are 12 variables 
in all. Using V tJ k and one is able to invert in order to 
obtain all the derivatives of the inverse metric in terms 
of the new variables: 



1 



(57) 



In addition to reducing the BSSN equations to first- 
order form, we also require the lapse function to be den- 
sitized, namely: 



-ae 6 * 



(58) 



where a is an arbitrarily prescribed function of space- 
time. The evolution equations (HJ, augmented by trivial 
evolution equations for Qi and V lJ k, read now: 



-2aA, 



(59a) 
(59b) 



Jirljsl V k,m - T^lkid ,j) +^7{7y i ,k 



7 



10 



30 



Qi = -—K,i + ... 



V ij k = -2a(A^, k ~8^A?>, a +U»n fk iA l ', a ^ 



(59c) 
(59d) 
(59e) 
(59f) 

(59g) 



The new evolution equations (|59t| and |o9g| are obtained 
simply by taking a time derivative of the definitions l|55a|) 
and i|55b|) and substituting the evolution equations (|8bl) 
and l)8a|) in the resulting right-hand sides. No linear com- 
binations with the constraints are used. 

In order to see that this initial value problem is 
strongly hyperbolic in the santadard sense , the eigen- 
value problem of the principal symbol in an arbitrary di- 
rection must be solved. This is equivalent to assuming 
that all the fields have a plane wave dependence of the 
form expi(£kX k — st) with real s and arbitrary £j so that 
£,i£,jl 13 — 1- The eigenvalue problem of the principal 
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To' 



(60a) 
(60b) 



(60d) 
(60e) 

(60f) 
(60g) 



The reader can verify that this problem has the follow- 
ing eigenvalues: 

• s = with 18 eigenvectors labeled by free val- 
ues of jij, 0,r\ the two components of Qi except 
J kl Qkt;i and the seven components of V 13 k other 
than 7 fc ^'fc6- 



s% = 





S(j) — 









sAij — 


-ote 40 (-7, 




7 

+ 30 7y " r 6 


sK = 








sf l = 






a 


sQi = 




v» k = 


-2a(i^ fe - 



• s = ±a with eigenvectors given by 

lij = 



r = 



4a 
"37 



A«6 



6s 



with free values of K and of the two components of 
A 13 other than A l °^j. These are six eigenvectors in 
all: three for s — a and three for s = —a. 



• s = ±ay 3/5 with eigenvectors given by 

lij 



r 

Qi 

K 

#3 



with free values of the component A kl ^k^i - T hese 
are two eigenvectors in all: one for s = cc-^/3/5 and 
one for s = —ay/3/5. 
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s = ±a^/7/10 with eigenvectors given by 



F 

Qi 
K 

A ij 



S \ 



s 

+ ~f j %lA lm tr, 

5 



with free values of the two components A* J £j other 
than A kl ^k£,i (which vanishes). These are four 
eigenvectors in all: two for s = ay7/10 and two 
for s = -ay/7/10. 

In all of the above, as in what follows, we are using the 
notation £ 4 = 7 y £j. 

All eigenvalues are thus real and there are 18 + 6 + 
2 + 4 = 30 linearly independent eigenvectors. The ini- 
tial value problem of Eqs. I|59|l is, thus, strongly hy- 
perbolic in the standard sense, and has six pairs of 
"non-static" characteristic fields traveling with non-zero 
speeds. Not all characteristic speeds are equal to the 
speed of light, though: there are only three pairs of 
incoming-outgoing characteristic fields traveling at the 
speed of light. The other three pairs of "non-static" char- 
acteristic fields travel at subluminal speeds. (This does 
not necessarily contradict , since in that reference the 
authors use a different first-order reduction of the BSSN 
equations.) In the following we write down explicitly 
the characteristic fields with non-vanishing characteris- 
tic speeds, with the aim of relating some of them to the 
boundary equations of Section ITTT1 

The three pairs of characteristic fields that travel at 
light speed s = ±a are 

± U K = K±&i k Q k (61) 

±U* = A ij - 2^A^ k H k + + i ij )A kl tei 

tk 



± 



i 



(62) 



Notice that f/'^Ty = ± E/ y £ J - =0. As observed previ- 
ously, they are labeled by free values of K and of the two 
components of A % i other than the projections along the 
characteristic direction £j. 

The pair of characteristic fields with s = ±ay/3/5 is 



4 f ^ 



v ml k ute k 



(63) 



As observed before, it is labeled by the projection 

Th e tw o pairs of characteristic fields with s — 
±a\/7/l0 are 



±jji 



(64) 



Notice that U' 1 ^ — 0. These are labeled by the two 
projections A % i£j other than A kl £k£h as anticipated pre- 
viously. 

In order to see which of these fields are prescribed by 
the boundary equations lJT5)l. (|T§1> . lt^|) and we re- 
strict attention to £j = (1,0,0)/ y/^ xx . Lowering the in- 
dex i with 7y we have 



z u z = 



A' 



At 



±lJ^ yk (f k -V xk ^ 



.r.r / ' \ 7 V 



.,xk , 



(65a) 
10 



2Vio 7 " fe v ~ i -^)^ A yY Qz 

(65b) 



Inverting for Ay and A x we have 



Al 



At 



-(-Uy + + Uy) 

-{-u z + +u z ) 



(66a) 
(66b) 



Clearly, the boundary equations H18[) and H19[) prescribe 
the time derivatives of the incoming fields _ U y and ~~ U z , 
respectively, in terms of their outgoing counterparts. 
Next, with & = (1, 0, 0)/Vt^, the fields ± U read 



j xx 3 V 3 



t rXX T 

- V k 



7 V x A sv kx 

--^=--e^Q k ^= 

15 y/j xx 3 \/~i xx 



which can be inverted as 



A xx Zj xx K 

3 1 



-(-U + +U) 



(67) 



(68) 



Furthermore, since A xx — ^j xx K = j x; > (A x — ^S X K) then 
we have: 



^(A X ~-5 X K) = ^(-U + +U) 
Isolating the term with j = x we have 



(69) 



r x {K - \k) = 'tl r u + +u) - r y A x y - r z A x z (?o) 
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Using H66fl in the right-hand side, we finally have 



%-\ K = e -Yru + + u)-^(-u y + ^u y ) 



;{-U z + +U z ) 



(71) 



Clearly, since the time derivatives of the incoming fields 
~~U y and ~U Z are already prescribed by (|T5j) and (JTDJ, 
then the boundary equation (12411 prescribes the time 
derivative of the incoming field ~U in terms of its outgo- 
ing counterpart. 

It only remains to understand the role of the fourth 
boundary equation l|2()l) . The reader can verify that the 
combination 
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~,xk/ 



(72) 



is one of the 18 characteristic fields with zero character- 
istic speed. But this is exactly equal to the combination 
that appears under the time-derivative sign in (|26|l up to 
terms of lower order, that is 

f x - 8T k Qk = 7 "((ln(det7 fei ))„ -T^Sm) + •• ■ (73) 

Clearly, thus, the fourth boundary equation l|26[) is a con- 
dition on the time derivative of the "static" characteristic 
field °U, and, as such, it should be identically satisfied. 
The fourth equation is not necessary in order to prescribe 
incoming fields. 

In summary, we have found that three of the six incom- 
ing characteristic fields are prescribed by three of the four 
equations G x a — 0. These are precisely the characteris- 
tic fields that travel at speeds different from the speed 
of light. The fourth equation is redundant, but must be 
satisfied by a static characteristic field as a consequence 
of the initial data and the evolution equations. 

This is entirely consistent with the analysis of Sec- 
tion llVI The results of this Section show that, in the case 
of the first-order reduction, there are only three boundary 
prescriptions in the set G x — 0, and the three incoming 
fields that travel at the speed of light remain freely speci- 
fiable. This is very similar to the case of the Einstein- 
Christoffel formulation discussed in and 0. The fields 
that remain freely specifiable are labelled by the compo- 
nents A v . and A^ — At, and the trace K of the extrinsic 

z y z 

curvature, as in the case of the Einstein-Christoffcl for- 
mulation. The exception is that in the present case it 
becomes clear that the incoming fields that remain free 
are physical in the sense that they are the only ones trav- 
eling at the speed of light (the same case cannot be made 
in the Einstein-Christoffel formulation because all its in- 
coming fields travel at light speed). 

It is not clear to us whether any of the conclusions 
regarding the incoming characteristic fields may be di- 
rectly transferred from the first-order reduction to the 
original BSSN equations, that is, Eqs. JSJ- The main use 
of the results of this Section is to add insight into the 
boundary-value problem of the (standard) BSSN equa- 
tions, as well as to increase our understanding of the role 



of the Einstein boundary conditions G x = in the initial- 
boundary- value problem of the Einstein equations across 
formulations. 



VI. CONCLUDING REMARKS 

We have written down explicitly the principal terms 
of the four components of the projection of the Einstein 
tensor perpendicularly to a timelike boundary in terms 
of the variables of the BSSN formulation, Eqs. ©■ The 
resulting expressions, Eqs. {TU, CHI, 12} and con- 
tain no second derivatives across the boundary of any of 
the fundamental variables, and no first derivatives across 
the boundary of any of the three first-order variables V . 
Consequently, setting them to zero results in four equa- 
tions internal to the boundary in the sense that they 
are constraints on the timelike boundary surface. Such 
boundary equations can then be interpreted as bound- 
ary conditions for the BSSN formulation of the Einstein 
equations. Trivially, these boundary conditions are nec- 
essary to the effect of satisfying the Einstein equations 
on the boundary surface. But far from trivially, we have 
demonstrated that three of them are necessary in order 
to ensure that the solution of the BSSN equations will 
satisfy the constraints in the future of the initial slice. 

A meaning to the three non-trivial boundary equations 
can be assigned as follows. A timelike boundary surface 
can be interpreted as an artificial boundary separating 
two regions in a spacetime generated by the same Cauchy 
surface (that is: a "big" spacetime with no boundary or 
with a null boundary). In this interpretation, the values 
on the timelike "boundary" surface are completely de- 
termined by evolution from the initial data. The initial 
data satisfy the momentum and Hamiltonian constraints. 
Therefore, the values on the timelike "boundary" surface 
must be related by some sort of evolution of the con- 
straints from their initial values. Three of the constraints 
propagate towards the "boundary" from data in the in- 
terior region of the Cauchy surface. But three of them 
propagate towards the "boundary" from Cauchy data in 
the other side of the timelike surface. These three "in- 
coming" constraints must affect the "boundary" values 
that are interpreted as incoming, and a prescription for 
their effect on the boundary values is needed. The point 
is that the three nontrivial equations in the set G x — 
accomplish such a prescription by realizing, in a sense, 
the evolution of the incoming initial constraints towards 
the "boundary" surface. 

Just like in the case of the ADM equations, the problem 
of the boundary conditions of the BSSN formulation re- 
mains open in the following sense. There are, in principle, 
15 variables whose boundary values are needed in order 
to proceed forward in an integration in time. Of these, 
four can be obtained from initial data by the four equa- 
tions G x — 0, which are, in fact, time-evolution equations 
internal to the boundary surface. It is not known, at this 
time, how many of the remaining 11 variables are freely 
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specifiable at the boundary, and how many are deter- 
mined by evolution from initial data. The potentially 
freely specifiable variables represent an indeterminacy in 
the sense that their role in shaping up the final solution 
is as important as that of the initial data. Thus, physical 
insight would be required to prescribe the (potentially 
free) boundary data in a manner analogous to the physi- 
cal insight required to determine the free initial data for a 
particular problem (e.g, for binary black hole collisions). 

We have also studied the related problem of a strongly- 
hyperbolic reduction of the BSSN equations, that is, 
Eqs. In this case, the boundary values that need 

to be specified are known to be those of six of the 30 
fundamental variables, of which three are prescribed by 
the Einstein boundary conditions. The others are deter- 
mined, in principle, from the initial values and the evo- 
lution equations. We have not concerned ourselves with 
the question of determining which linear combinations of 
G x a = lead to a well-posed initial-boundary-value prob- 
lem, since our present first-order reduction is meant only 
as an auxiliary device in order to shed light on the orig- 
inal problem and is not expected to be used for numeri- 
cal simulations. The problem of separating the Einstein 
boundary conditions for this first-order reduction into a 
well-posed set and an ill-posed set remains open at this 
time. 



Other issues that remain wide open are whether or 
not the (standard) BSSN formulation with the Ein- 
stein boundary conditions can be implemented in a nu- 
merically stable manner in the rigorous sense of finite- 
differences, and whether the Einstein boundary condi- 
tions have any effect on the run time of numerical sim- 
ulations. Because the BSSN formulation is not strongly 
hyperbolic in the standard sense, the former issue lies 
outside of our current interests. The latter issue is also 
of fundamental importance to numerical simulations be- 
cause it has been observed that constraint violations play 
a crucial role in shortening the run time 16]. Our ex- 
pectation is that the Einstein boundary conditions, be- 
ing "constraint-preserving" in the sense discussed in Sec- 
tion can potentially increase the run time of numeri- 
cal simulations, assuming that other potential numerical 
instabilities are under control. 



Acknowledgments 

This work was supported by the NSF under grants No. 
PHY-0244752 to Duquesne University, and No. PHY- 
0135390 to Carnegie Mellon University. 



[1] S. Frittelli and R. Gomez, Class. Quantum Grav. 20, 
2379 (2003), gr-qc/0302032. 

[2] S. Frittelli and R. Gomez, Phys. Rev. D 68, 044014 
(2003), gr-qc/0302071. 

[3] S. Frittelli and R. Gomez, Einstein boundary con- 
ditions in relation to constraint propagation for the 
initial-boundary value problem of the Einstein equations, 
gr-qc/0310064. 

[4] A. Anderson and J. W. York, Phys. Rev. Lett. 82, 4384 
(1999). 

[5] J. W. York, in Sources of Gravitational Radiation, edited 
by L. Smarr (Cambridge University Press, Cambridge, 
1979). 

[6] T. W. Baumgarte and S. L. Shapiro, Phys. Rev. D 59, 
024007 (1999). 

[7] M. Shibata and T. Nakamura, Phys. Rev. D 52, 5428 

(1995) . 

[8] S. Frittelli and O. A. Reula, Phys. Rev. Lett. 76, 4667 

(1996) . 

[9] U. Sperhake, K. L. Smith, B. Kelly, P. Laguna, and 



D. Shoemaker, Impact of densitized lapse slicings on evo- 
lutions of a wobbling black hole, gr-qc/0307015. 

[10] B. Bruegmann, W. Tichy, and N. Jansen, Numerical sim- 
ulation of orbiting black holes, gr-qc/0312112. 

[11] P. Laguna and D. Showmaker, Class. Quantum Grav. 19, 
3679 (2002). 

[12] B. Gustaffson, H.-O. Kreiss, and J. Oliger, Time- 
dependent problems and difference methods (Wiley, New 
York, 1995). 

[13] G. Nagy, O. Ortiz, and O. A. Reula, Strongly hy- 
perbolic second order Einstein's evolution equations, 
gr-qc/0402123. 

[14] S. Frittelli and O. A. Reula, J. Math. Phys. 40, 5143 
(1999). 

[15] O. Sarbach, G. Calabrese, J. Pullin, and M. Tiglio, Phys. 
Rev. D 66, 064002 (2002). 

[16] L. E. Kidder, M. A. Scheel, S. A. Teukolsky, E. D. Carl- 
son, and G. B. Cook, Phys. Rev. D 62, 084032 (2000). 



